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Jeff Greensite 



1. Introduction 



From time to time there have been efforts [1-13] to obtain the Yang-Mills vacuum wavefunc- 
tional, often in lower dimensions, to see if anything can be learned about confinement and the mass 
spectrum. The problem is to solve for the ground state To [A] of the Yang-Mills Hamiltonian H, 
which has its simplest form in temporal gauge 
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The price for this simplicity, in temporal gauge, is that all physical states are required to satisfy a 
constraint 
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which enforces invariance under infinitesimal gauge transformations. 

Our claim is that the ground state solution in D = 2 + 1 dimensions, in temporal gauge, is 
approximated by 1 



T [A]=exp 
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(1.3) 



where B" = F" 2 is the color magnetic field strength, D 2 is the co variant Laplacian in the adjoint 
color representation, Xq is the lowest eigenvalue of — D 2 , and m is a constant of order g 2 . In support 



of this claim, we will argue that (1.3) (i) is a solution of the YM Schrodinger equation in the g — > 
limit; (ii) solves the YM Schrodinger equation in the strong field, zero-mode limit; (iii) confines 
if m > 0, and that m > seems energetically preferred; and (iv) results in the numerically correct 
relationship between the mass gap and the string tension. A more complete report of our results 
can be found in a recent article [15]. 



2. Free-Field and Zero-Mode Limits 



In the free-field g = limit, eq. (1.3) becomes 



T [A]=exp 



\ J d 2 xd 2 y(d l A a 2 {x) - d 2 Al{ 



d l A h 2 {y)-d 2 A\{y) 



(2.1) 



which is the known ground state solution in the abelian, free-field case, and an obvious starting 
point for any investigation such as ours. Apart from the free-field limit, there is another, quite 
different limit which can be treated analytically. Consider gauge fields which are constant in space, 
but variable in time, in D=2+l dimensions. The Lagrangian is 



d 2 x 



d,A k ■ d,A k - g 2 (A l x A 2 ) • (Aj x A 2 ) 



V 



d t A k ■ d t A k - g 2 (A i x A 2 ) • (A i x A 2 ) 



(2.2) 



'Closely related proposals have been made in the past by Samuel [7] and Diakonov [14], cf. the discussion in 
ref. [15]. 
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with Hamiltonian operator 

H 



4^ + ^ V(AlXA2) - (AlXA2) 



(2.3) 



where V is the volume of 2-space. With some algebra, one may verify that 

(Ai xA 2 )-(Ai xA 2 ) 



= exp 



VN 2 + |A 2 | 2 



(2.4) 



satisfies the zero-mode Yang-Mills Schrodinger equation up to 1 /V corrections. For comparison, 
we consider our proposed vacuum state ( |L3| ) in the strong A-field limit, where the covariant Lapla- 
cian is dominated by the gauge-field zero-mode, i.e. 



-D 2 )% = g z 8(x-y) (A 2 +A 2 )8 ah -A\A\-A^ 2 



(2.5) 



In this limit, it is not hard to show that *^o[A] of eq. (1.3) precisely reduces to (2.4). Thus the 
proposed ground state not only has the correct free-field limit, but also agrees with the calculable 
ground state of the zero-mode Schrodinger equation in an appropriate strong-field limit. 



3. Dimensional Reduction and Confinement 



A long time ago it was suggested [1] that at large distance scales, the pure Yang-Mills vacuum 
in a confining theory looks like 



eff , 



exp 



(3.1) 



This vacuum state has the property of dimensional reduction: Computation of a spacelike loop 
in d + 1 dimensions reduces to the calculation of a Wilson loop in Yang-Mills theory in d Eu- 
clidean dimensions. In d = 2 dimensions the Wilson loop can be calculated analytically, and we 
know there is an area-law falloff, with Casimir scaling of the string tensions. Now suppose, in 



the proposed vacuum (1.3), we expand the B-field in eigenmodes of the covariant Laplacian, i.e. 
B a {x) = Y*nbn§n{ x ) where — D 2 <p t "(x) = X„(j)%(x). Define the "slow" component 



(3.2) 



n=0 



with mode cutoff n max defined such that A„ mar — Ao <C m 2 . Then the portion of the squared wave- 
functional ( |0| ) which is quadratic in B slow , to leading order in l/m 2 , is 



exp 



- / d 2 xB slow B skm 
m 



(3-3) 



which has the dimensional reduction form. The string tension for fundamental representation Wil- 
son loops is easily computed in two Euclidean dimensions, and in lattice units it is a = 3m/ (4/3). 
Suppose we turn this around, and choose the mass parameter m = |j3a. Then our proposed vacuum 
wavefunctional must imply a definite mass gap, which we would like to calculate. 
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4. The Mass Gap 

To get the mass gap, we need to compute the connected correlator 

(B 2 (x)B 2 (y)) conn = ((B a B«) x (B b B h ) y ) - ((B a B a ) x ) 2 (4.1) 
in the probability distribution 



P[A} = \y [A}\ 2 = exp 

where 

rah 
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ab 



K%\A] = - (4.3) 

Numerically this looks hopeless! Not only is the kernel K 1 ™ highly non-local, but it is not even 
known explicitly for arbitrary gauge fields. We have, nonetheless, found a way of carrying out the 
numerical simulation, which relies on the fact that, after gauge-fixing, the variance in the kernel 
among thermalized configurations is negligible. The method is described in ref. [15]; here we will 
only present the results. 

Observables of interest include the eigenvalue spectrum {A,,} of the adjoint covariant Lapla- 
cian, and the connected field-strength correlator (B 2 (x)B 2 (y)) conn <x G(x — y), where 

G(x-y) = ({K-XXK- l )t) (4.4) 



Choosing m = gjSa, the mass gap can be extracted from G(R). For comparison, we can also 
compute these observables on time-slices of lattices generated by ordinary lattice Monte Carlo in 
D = 3 Euclidean dimensions. We will refer to these time-slices as MC lattices, and they can be 
regarded as drawn from the probability distribution |*Pq [U}\ 2 , where *Pq is the ground state of the 
transfer matrix of the Wilson action. Lattices generated by the method described in ref. [15], which 
simulates our proposed ground state wavefunctional, will be referred to as recursion lattices. 

Fig. l(a)| is a plot of eigenvalue vs. mode number for the operator — D 2 — Xq + m 2 , at j8 = 18, 



from ten independent 50 x 50 recursion lattices (no averaging). Also plotted, but indistinguishable 
from the other spectra, is the rescaled spectrum of the large- volume zero-field operator —V 2 + m 2 . 
It is obvious from this figure that there is very little variance in the spectrum of — D 2 — Ao from one 
thermalized lattice to the next, moreover, this spectrum is very close to that of a free-field theory. 



Fig. |l(b)| shows our data for G(R), obtained from ten recursion lattices, and ten MC lattices. Note 
the tiny values of G(R) ~ 10~ 12 obtained at the larger R values. This requires a near-absence of 
fluctuation in K from one thermalized lattice to the next. Note also that G(R) obtained on MC 
and recursion lattices agree very closely with one another. The mass gap M is extracted from a 
two-parameter (c,M) fit of the recursion lattice data to the function 2 

-MR 

G fit (R)=c(\ + \MR) 2 ^- (4.5) 

The best fit to the data for G(R) at j3 = 18 is shown in Fig. |2(a)| . Our final results for the mass gaps 
at a variety of j8 values are shown in Fig. 2(b)| , together with the results for the + glueball obtained 



by Meyer and Teper [16], using standard methods. The agreement is clearly very good. 
2 This expression is motivated from the functional form of K^zK^ where Kq 1 = \/— V 2 +m 2 . 
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Figure 1: Simulations of the vacuum state at j3 = 18 and lattice extension L = 50. (a) Eigenvalue spectra of 
the operator — D 2 — Ao + m 2 . Also plotted is the spectrum of the large-volume zero-field operator — V 2 +m 2 . 
(b) The correlator G(R) computed on recursion and MC lattices. 
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Figure 2: (a) Best fit (dashed line) of the recursion lattice data for G(R) by the analytic form given in 
eq. (4.5). (b) Mass gaps extracted from recursion lattices at various lattice couplings, compared to the + 
glueball masses in 2+1 dimensions obtained in ref. [16] (denoted "expt") via standard lattice Monte Carlo 
methods. Errorbars are smaller than the symbol sizes. 



5. Why confinement? 

We have seen that there is a mass gap if K 1 has a finite range, and a non-zero string tension 
(by the dimensional reduction argument) if K is also finite range. Both of these conditions are 
obtained if the mass parameter m is non-zero, so the question "why confinement?" boils down to 
the question: why is m non-zero? An obvious approach is to treat m as a variational parameter, and 
use it to minimize (H) . To simplify matters, noting the negligible fluctuation in K, we make the 
drastic approximations of (i) neglecting functional derivatives of the kernel K in computing (H); 
and (ii) ignoring correlations between B and K. With these approximations, the VEV of H turns 
out to be 

(H) = \ /W-D2-Ao + m2 + ^Tr tf£ =) (5.1) 
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R[U] = 




contours 



I larger contours | l L !L » 1-========- ; 

(a) (b) 

Figure 3: (a) The first few terms in the strong-coupling expansion of the lattice vacuum state M'oP^'L with 
R[U] = logroll/]), (b) How 1x2 rectangles in R[U] screen an adjoint Wilson loop. The adjoint Wilson 
loop is denoted by a heavy solid line. The overlapping 1x2 rectangles are indicated by (alternately) light 
solid and light dashed lines. 



There is a competition between the two terms on the rhs; the first increases with m 2 , while the 
latter decreases. In an abelian theory, where D 2 = V 2 and Ao = 0, it is trivial to show that (H) is 
minimized at m = 0, implying no confinement, and no mass gap in the abelian case. For a non- 
abelian theory matters are different, chiefly because Ao > 0. A naive calculation, which ignores 
any dependence of the {A„} on m 2 , finds that m 2 = (Ao) at the minimum. While a more accurate 
calculation would no doubt change somewhat the numerical value of m 2 at the minimum, it seems 
unlikely that this value would be shifted exactly to zero. 



6. What about N-ality? 



Dimensional reduction at large scales, from three to two dimensions, implies Casimir scaling 
of the string tensions; in particular, the adjoint string tension is non-zero. While this is correct at 
intermediate distance scales, asymptotically the string tension should depend only on the N-ality 
of the representation. The transition from Casimir scaling to N-ality dependence must be coming 
from corrections to dimensional reduction. 

It is worth examining how N-ality dependence is achieved when the vacuum state is com- 
puted by the strong-coupling method developed in ref. [2]. In this case the ground state has the 
form *^o[^] = exp [/?[[/]], where R[U] is a sum over contours (plaquette, 1x2 rectangle, etc.) con- 
structed from p plaquettes, with each contour multiplied by a constant of order (5 2p . The leading 



contributions, at strong-coupling, are shown in fig. p(a)| . It is not hard to show that the perimeter 
law for large adjoint loops is obtained, in the strong-coupling expansion of ^[t/], by tiling the 



perimeter of the loop with overlapping 1x2 rectangles, as shown in Fig. |3(b) . 

It is interesting that the rectangle term not only screens adjoint loops, but is also responsible 
for the leading correction to dimensional reduction. An expansion of the strong-coupling vacuum 
in powers of lattice spacing yields [17] 



¥ [t/]=exp 



d 2 x {qKqB 2 -a i K 2 B{-D 2 )B+ ...) 



(6.1) 
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The first term represents dimensional reduction, and all four contours shown in Fig. 3(a)| contribute 



to Kq. Only the rectangle contour, however, contributes to the second term, which supplies the 
leading correction to dimensional reduction. It is also the rectangle contour which couples together 
field strengths in neighboring plaquettes. 

Returning to our proposed vacuum state and expanding the kernel in powers of 1/m 2 , the part 
of *Pq quadratic in B sIow , up to next- to-leading order in 1/m 2 , is 



exp 



- [d 2 x ( B sl0W B sl0W - B slow ° 2 = ^ B slow + ... 

2m 1 



m 



(6.2) 



Note that the correction to the first (dimensional reduction) term is similar to, and has the same sign 
as, the leading correction to dimensional reduction in the strong-coupling vacuum state, generated 
by the rectangle term. Our speculation is that, as in strong coupling, this second term is associated 
with the screening of higher-representation color charges, converting Casimir scaling to N-ality 
dependence of the higher-representation string tensions. 
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